■5^-y 7 

p. 7 

Remote Sensing Estimates of Cirrus Particle Size for Tropical 
and Midlatitude Cirrus : Hexagonal Crystals 

Bruce A. Wlellcki and Patrick Mlnnis 
Atmospheric Sciences Division 
NASA Langley Research Center 
Hampton, VA 23681-0001 

Robert Ardulnl and Lindsay Parker 
Lockheed Engineering and Sciences Corp. 

Hampton. VA 23666 

Sl-Chee Tsay 

Universities Space Research Association 
NASA Goddard Space Flight Center 
Greenbelt MD 20771 

Yoshlhlde Takano and Kuo-Nan Llou 
Department of Meteorology, University of Utah 
Salt Lake City, Utah 


and Ice Spheres. 

V 

N94-22342 


1. INTRODUCTION 

A large discrepancy exists in current estimates of a mean cirrus particle size 
appropriate for calculations of the effects of these ice clouds on solar and thermal infrared 
radiative fluxes. For spheres with large size parameter (x = 2xr/A,> 30, where r Is particle 
radius), and moderate absorption (n* x < 1 , where n* Is the imaginary index of refraction for 
Ice), the optimal effective particle radius (Hansen and Travis, 1974) Is given by: 

r, = jr 3 n(r)dr/jr 2 n(r)dr. 

For the remote sensing of cirrus particle size at wavelengths of 0.83, 1.65, and 2.21pm, a 
50pm ice sphere would have a size parameter of about 200, and values of n 1 x of 0, 0.045, and 
0.06, satisfying the above conditions. 

However, while r e Is a well-defined parameter for spheres, this cross-section area- 
weighted particle radius can only be extended to non-spherical particles by defining some 
equivalent sphere, typically an equivalent volume or equivalent cross-section area sphere. 
Using equivalent volume spheres, values of r e obtained over Lake Michigan on October 28, 
1986, during FIRE phase I varied from 200pm (King Air 2D Imaging probes) to 60pm (Landsat 
reflectances at 0.83, 1.65, and 2.2pm, Wlellcki et al., 1990), to 25pm (HIS spectrometer thermal 
emission between 8 and 12pm, Ackerman et al., 1990). Three major uncertainties were 
identified in this comparison: small ice particles missed by the 2D-C aircraft probes, uncertain 
ice refractive index, and uncertainties in the single scatter albedos and scattering phase 
functions used In the radiative calculations. 

Since the first FIRE cirrus results, advances have been made In all three areas. The 
present paper reports on improvements In the radiative modeling of ice particles at 0.83, 1.65 
and 2.2 1 pm wavelengths appropriate for comparisons to Landsat Thematic Mapper data. 

The paper also Includes new results for Landsat observations of ice clouds in the eastern and 
western tropical Pacific. 

2. RADIATIVE MODEL 

Radiative transfer calculations used to predict Landsat spectral radiances are 
performed as outlined in Wlellcki et al., 1990 with two exceptions. First, the finite difference 
radiative transfer method has been replaced by a discrete ordinate method, hereafter termed 
DISORT (Stamnes et al., 1988; Tsay et al., 1991). The overall accuracy of these two methods. 
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however, are similar. Second, a hexagonal Ice crystal model (Takano and LIou, 1989) Is used 
in the present study for both the single scatter albedo (o>), and scattering phase function (<I>) 
used in the DISORT calculations. In the earlier study by Wicltckl cl al.. (1990) Mle spheres 
were used to determine 10 , while <I> was taken from Jabora lory measurements of ice crystals by 
Volkovitskiy et al. 1980, hereafter denoted as VPP. The VPP phase function was only available 
for X = 0.63pm, but was used for all three channels, in the absence of better information. The 
present study uses phase functions derived for randomly oriented hexagons of width 2a and 
length L (see Table 1) determined using ray tracing at X = 0.55, 1.6, and 2.2pm (Takano and 
Llou, 1989). Figure 1 shows the hexagonal phase functions for L/2a = 120/60 at all three 
wavelengths. The figure demonstrates the importance of absorption in reducing the amount 
of side scattered radiation, as the average imaginary Index of refraction for ice increases with 
wavelength. Table 1 summarizes the radiative modeling differences between the two studies. 


TABLE 1. 


Radiative Modeling Changes 


Radiative Model Scattering Phase Function 

Current Results: 

Discrete Ordinate Hexagons, ray tracing: 2a/L = 

20/20, 50/40, 120/60, 
300/100, 750/160 pm/pm; 

X = 0.55, 1.6, 2.2 pm 

Wlellcki et aL. 1990: 

Finite Difference Laboratory: 20 pm plates: 

X = 0.6 pm only 


Single Scatter Albedo 

Hexagons, ray tracing: 2a/L = 
20/20, 50/40, 120/60, 

300/100. 750/160 pm/pm; 

0.83, 1.56 - 1.8, 2.10 - 2.35 pm 

Mie : r e = 4, 8, 16, 32, 64, 128pm 
0.83, 1.56 - 1.8, 2.10 - 2.35 pm 


Figure 2 compares the new hexagonal crystal results for the Landsat 0.83 pm and 1.65 pm channels 
with the previous VPP/Mie sphere results at a solar zenith angle of 45°. Calculations were performed 
at 0.83pm extinction optical depths of 0. 1, 0.2, 0.5, 1, 2, 3.5, 5, 10, 20, and 50 and are marked by 
symbols In Figure 2. Takano and Liou (1989) and others have shown that particle absorption 
efficiency Qabs scales with particle volume. In order to compare with the earlier Mie sphere estimates 
of Qabs< (which for large particles with 2nr/X » 1 determines to) we have given the radius of an 
equivalent volume sphere for each hexagonal crystal size. Figure 2 shows that for large optical 
depths the new calculations give results similar to an equivalent volume sphere. At smaller optical 
depths, however, the changes in the hexagonal scattering phase function with X (Figure 1) cause 
larger deviations from the earlier results, so that the hexagonal crystal rvol is as much as a factor of 
2 smaller than the VPP/Mie sphere results. In summary, at large optical depths the single scatter 
albedo dominates the result, and equivalent volume spheres give reasonable results, while at optical 
depths less than about 10, the X dependence of the scattering phase function plays a large role. 


3. RESULTS 

Figure 3 compares the hexagonal crystal predictions for Landsat reflectances at 
0.83|im and 1.65pm with Landsat observations. Landsat observations are shown at full 
spatial resolution of 30 meters, and are uniformly sampled over a set of independent Landsat 
cloud scenes. Each Landsat scene is a 60 km by 60 km region. The observations are divided 
into three sets by geographic regions. Figure 3a gives results for 14 Landsat (60 km) 2 scenes 
taken from the western tropical Pacific near the TOGA/COARE region (13S - 13N, 150E-175E). 
Figure 3b gives results for 6 scenes from the Eastern Pacific (ON - 13N. 95W-135W). Figure 3c 
gives results for 5 scenes for northern mid latitudes (30N - 60N, 60W-130W). 

Figure 3 shows several interesting characteristics. For all 25 regions, greater than 
98% of the observations fall between the 50/40 and 300/100 hexagonal crystal results (rvol 
between 23 and 78 pm). The median value agrees well with the 120/60 crystal size (r V ol = 41 
pm). The Eastern Pacific results appear to be bimodal (but no single scene shows both 
modes). The mid latitude cases and most of the optically thin tro pi cal cirrus (0.83 pm 
reflectance below 0.2) show a smaller size range of about 50/40 to 120/60. For all data, 
however, very few observations are found with particle sizes smaller than 50/ 40. In fact, the 
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small number of observations following the 20/20 results for the mid latitude cirrus in Figure 
3b are supercooled water droplets from the 10/28/86 case in WIelicki et al. (1990). 

4. CAVEATS 

The great complexity of ice crystal shapes in clouds indicates that even hexagonal 
crystal interpretations of remotely sensed particle size must be viewed with caution. In 
particular, studies of the dependence of the scattering phase function on particle shape and 
imaginary index of refraction are especially important. Final resolution of these Issues will 
require consistency of aircraft, and various remotely sensed particle estimates including near- 
infrared and thermal Infrared radiance techniques, as well as new techniques using multiple 
wavelength active Ildar and radar systems. It should also be noted that the particle sizes 
Inferred from Figure 3 would represent a vertical average over the entire cloud depth for 
optically thin cloud, while the sizes would represent the microphysical conditions nearer 
cloud top (optical depth 3-6) for optically thick cloud. 

Even with these difficulties, the current results suggest that there might be strong 
similarities in mid latitude and tropical Ice microphysics, at least for optically thin cirrus. 
Additional cases of mid latitude ice clouds with larger optical depths are needed. 
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Figure 3a 
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Figure 3b 




